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DYNAMIC STABILITY DERIVATIVES OF A TWIN- JET FIGHTER MODEL 
FOR ANGLES OF ATTACK FROM -10° TO 110° 

By Sue B. Grafton and Charles E. Libbey 
Langley Research Center 

SUMMARY 

A low-speed investigation was conducted in the Langley full-scale tunnel to deter- 
mine the dynamic stability derivatives in pitch, roll, and yaw over an angle-of-attack 
range of -10° to 110° for a twin- jet swept-wing fighter airplane. The study consisted of 
forced- oscillation tests of a 0.13- scale model of the airplane at a Reynolds number of 
1.33 X 10®. Several oscillatory frequencies and amplitudes were investigated to deter- 
mine the effects of these variables on the stability derivatives. The effects of the verti- 
cal and horizontal tail, and of the horizontal -tail deflection, on the derivatives were also 
evaluated. 

The results of the investigation indicated that the model exhibited stable values of 
damping in pitch over the entire angle-of-attack range, but marked reductions of damping 
in roll were measured at the stall, and unstable values of damping in yaw were present 
for the very high angles of attack associated with flat spins. These unstable values of 
damping in yaw, which constitute propelling moments in a spin, were produced by aero- 
dynamic interference between the vertical and horizontal tails. Either removal of the 
horizontal or vertical tail or full up deflection of the horizontal tail eliminated the pro- 
pelling moments. 


INTRODUCTION 

The National Aeronautics and Space Administration is currently engaged in a 
research program to develop and validate theoretical methods for prediction of airplane 
stall and spin characteristics. A major portion of this program involves correlation 
between theoretical calculations and data obtained by free-flight tests using dynamically 
scaled radio-controlled models of fighter airplanes. Previous theoretical and experimen- 
tal results obtained for a variable -sweep fighter with a long pointed nose are reported in 
reference 1. 

The present investigation was conducted to measure the dynamic stability deriva- 
tives of a twin-jet swept-wing fighter model over the angle-of-attack range associated 
with spinning; these data are intended to serve as aerodynamic input for additional 


theoretical studies of the stall and spin characteristics of this particular configuration. 
The investigation consisted of forced- oscillation tests which were conducted over an 
angle -of -attack range from -10° to 110° and included the effects of frequency and ampli- 
tude of the oscillatory motion. Tests were also conducted with the vertical and horizon- 
tal tails removed. The effect of horizontal-tail deflection on the yawing dynamic deriva- 
tives is also presented. Pertinent static force tests were conducted to aid in the analysis 
of the dynamic data. 

SYMBOLS 
Physical Concepts 

The longitudinal and lateral-directional data presented herein are referred to the 
body-axis system (see fig. 1). All data are referred to a center -of -gravity position of 
33-percent mean aerodynamic chord. 

In order to facilitate international usage of data presented, dimensional quantities 
are presented both in the U.S. Customary Units and in the International System of Units 
(SI). Equivalent dimensions were determined by using the conversion factors given in 
reference 2. 
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My 


pitching moment, ft- lb (m-N) 


M 2 yawing moment, ft- lb (m-N) 

p rolling velocity, rad/sec 

q pitching velocity, rad/sec 

free-stream dynamic pressxire, Ib/ft^ (N/m^) 

r yawing velocity, rad/sec 

5 wing area, ft^ (m2) 

u,v,w components of resultant velocity along X, Y, and Z body axes, 

respectively, ft/sec (m/sec) 

V free-stream velocity, ft/sec (m/sec) 

resultant linear velocity, ft/sec (m/sec) 

X,Y,Z body reference axes 

a angle of attack, deg or rad 

jS angle of sideslip, deg or rad 

6j^ horizontal -tail deflection, positive when trailing edge down, deg 

p air density, slugs/ft3 (kg/m3) 

6 angle of pitch, deg or rad 

A0 pitch- angle increment, deg 

(j) angle of roll, deg or rad 

A(p roll-angle increment, deg 
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angle of yaw, deg or rad 


A\{/ yaw-angle increment, deg 

CO angular velocity, 2jrf, rad/sec 

Dot over symbol indicates derivative with respect to time. 

Coefficients and Derivatives 

Results presented herein are given in terms of the coefficients and derivatives 
defined in the tabulations that follow: 
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In the present investigation the term "in-phase derivative” refers to any one of the 
oscillatory derivatives that is based on the components of forces and moments in phase 
with the angle of pitch, roll, or yaw produced in the oscillatory tests. The term "out-of- 
phase derivative" refers to any one of the oscillatory derivatives that is based on the 
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components of forces and moments 90° out of phase with the angle of pitch, roll, or yaw. 
The oscillatory derivatives of the present investigation were measured in the following 
combinations : 



Pitching 


Rolling 


Y awing 


Out of phase 

Cmq + ^m^j, 

CfJ + C-NT . 

iNq iNq, 

C A ^ A * 

^Yp + Cy^ sin 0! 
Cnp + sin a 

C^p + sin Of 

CYr " ^Y^ cos a 
Cur - cos a 

CZr - « 




APPARATUS AND MODEL 
Wind Tunnel 

The tests were made in the 30- by 60-foot (9.1- by 18.3-m) open-throat test sec- 
tion of the Langley full-scale tunnel with the model mounted about 10 feet (3.05 m) above 
the groundboard with its wings in a vertical plane. The model was so small in proportion 
to the tunnel test section that no wind-tunnel jet-boundary or blockage corrections were 
required. 
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Model 


The investigation was conducted with a 0.13-scale model of a two- place twin-jet 
high-performance fighter designed for land and carrier-based operations. A three- view 
sketch showing the general layout of the model is presented in figure 2; a photograph of 
the model is shown in figure 3, and geometric characteristics of the full-scale airplane 
are listed in table I. The longitudinal-control system of the configuration consists of an 
all-movable horizontal tail which incorporates negative dihedral of 23° to satisfy the lon- 
gitudinal stability requirements in the normal operational flight range . Lateral control 
is provided by spoilers as well as ailerons. The ailerons deflect downward while the 
spoilers deflect upward. The left aileron and right spoiler operate simultaneously as do 
the right aileron and left spoiler. The directional-control system consists of a conven- 
tional rudder. The maximum control-surface deflections are as follows: 


Rudder deflection, deg ±30 

Horizontal-tail deflection (trailing edge), deg 21 up, 9 down 

Aileron deflection, deg 0 up, 30 down 

Spoiler deflection, deg 45 up, 0 down 


All dynamic force tests were made with a single-strut support system and an inter- 
nal six-component strain-gage balance. The test setups for pitching, rolling, and yawing 
are illustrated in figure 4 and the equipment and readout system is described in refer- 
ence 3. The model and the strain- gage balance were mounted on the oscillation sting 
assembly so that the moment reference center of the balance was at the center-of-gravity 
location shown in figure 2 (33 percent c) and was on the axis of rotation for all test condi- 
tions. Oscillatory motion was imparted to the model by means of a flywheel that was 
driven by a 3-horsepower (2.2-kilowatt) variable-speed electric motor and a system of 
pushrods and bellcranks . The amplitude of the oscillatory motion was adjusted by varying 
the location of the lower pivot point of the vertical connecting rod along the radius of the 
flywheel. The frequency of the oscillatory motion (limited to about 2 cycles per second) 
was varied by changing the speed of the electric motor. 

A precision sine-cosine potentiometer, which generated voltage signals proportional 
to the sine and cosine of the flywheel rotation angle, was coupled directly to the flywheel 
shaft and provided electrical signals proportional to the angular displacement of the model. 
These signals were used in the data readout procedure described in detail in reference 3. 

TESTS 

The forced-oscillation tests in pitch, roll, and yaw were conducted over an angle-of- 
attack range from -10° to 110°. The tests were conducted in the Langley full-scale tunnel 
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at a speed of 100 feet per second (30.48 m/sec) which corresponds to a Reynolds number 
of about 1.33 X 10® based on the mean aerodynamic chord of the wing. Tests were made 
for the complete model and also for the model with the horizontal or vertical tails indi- 
vidually removed, and with the horizontal tails deflected to -21°. The tests were con- 
ducted for amplitudes ranging from ±2.47° to ±10.5° and at frequencies of 0.5, 0.7, 1.0, 
and 1.3 and 1.5 cycles per second over the entire angle -of -attack range. 

RESULTS AND DISCUSSION 

The more pertinent results with regard to stall and spin characteristics have been 
reviewed in previous reports. Lateral-directional characteristics of this airplane at the 
stall have been reported in reference 4 while aerodynamic factors affecting flat spin ten- 
dencies are presented in reference 5. The results presented herein, which form a com- 
plete documentation of dynamic stability derivatives at low speed, are therefore presented 
with a minimum amount of analysis. 

Static Aerodynamic Characteristics 

Longitudinal c haracteristics.- The static longitudinal characteristics of the model 
as functions of angle of attack are presented in figure 5. Results are shown for 
horizontal-tail deflection angles of 0° and -21° (corresponding to stick full back), and 
also for the horizontal-tail- off condition. The normal-force coefficients presented in 
figure 5(a) indicate that significant flow separation and stall began to appear at an angle 
of attack of about 12°. The variation of the axial-force coefficient shown in figure 5(b) 
also indicates onset of flow separation at about 12°. The results presented in figure 5(c) 
indicate that the model is statically unstable over a small angle -of -attack range just at 
the onset of the stall (12° to 18°), but has control effectiveness sufficient to overcome 
any pitchup or deep-stall trim condition. With the stick full back = -21°), the model 
can be trimmed at an angle of attack well above the stall, and the model is statically stable 
at all elevator positions throughout the spin angle -of -attack region. Above the stall, 
pitch-control effectiveness decreases to a minimum at about a - 60° but then increases 
as angle of attack approaches 90°. The results obtained with the horizontal tail off indi- 
cate that the horizontal tails contribute a very large diving moment in the spin angle-of- 
attack region. These results also reveal that the small region of static instability near 
the stall is also present with the horizontal tail off and is apparently related to wing stall. 

Lateral- directional characteristics.- The static lateral-directional stability deriva- 
tives of the model are presented in figure 6 as functions of angle of attack based on values 
of the coefficients at jS = ±5°. Results are shown for horizontal-tail deflections of 0°, 
-21°, and (over a limited angle-of-attack range) for the horizontal-tail-off condition. 

These results indicate that the model is directionally stable up to angles of attack slightly 


8 



above the stall (about 21°) and that effective dihedral increases about linearly with angle 
of attack up to a = 12°, where flow separation begins, and causes a marked reduction in 
Cl until it becomes zero at about the same angle of attack at which the model becomes 

directionally unstable. This combination of directional instability and zero effective 
dihedral at angle of attack in the region just beyond the stall is conducive to directional 
divergence problems which might lead to an inadvertent spin entry. Reference 4 dis- 
cusses the factors which cause the directional divergence to occur for this airplane con- 
figuration. The results of reference 4 indicated that loss of directional stability resulted 
from a combination of an adverse sidewash region at the rear of the model and a reduc- 
tion in dynamic pressure at the vertical-tail location. 


For the angle -of -attack range above 40° the model with controls neutral ^6^ = 0°j 
remained directionally unstable. Deflection of the horizontal tail to the full-back-stick 
condition (Sh = -21°) caused the directional instability above the stall to be more severe 
up to angles of attack of about 60°, and again at angles of attack above 88°. This result 
indicates the probability of an interference effect between the horizontal and vertical tails 
at these angles of attack. The model had positive effective dihedral for the entire angle- 
of-attack range above a = 40° and was not appreciably affected by deflection or 

removal of the horizontal tail in the angle -of -attack range. A discussion of the signifi- 
cance of the aerodynamic characteristics at very high angles of attack in terms of the flat 
spin is given in reference 5. 


Dynamic Stability Derivatives 

Pitching.- Comparisons of the effects of amplitude and frequency on the out-of- 
phase and in-phase pitching derivatives are presented in figures 7 and 8, respectively, 
for a limited angle -of -attack range (-7° to 28°). The out-of-phase pitching derivatives 
for the two largest amplitudes tested are presented in figure 9 for the entire angle-of- 
attack range investigated (-12° to 108°). 

The results presented in figures 7 and 8 indicate that damping-in-pitch increases 
in the a-range where static instability occurred (12.5° to 17.5°) in that the derivative 
combination Cjnq + became more negative. This result is probably related to 

flow separation on the wing. The effect of removal of the horizontal tails was to reduce 
Cmq + Cm^ in the low angle -of -attack range, but this increment was gradually reduced 
above the stall. The in-phase derivative - k2Cmq) reflects the trend of the static 

data; that is, the model exhibited reduced stability for 12° to 18° and this reduction was 
caused by the wing-fuselage combination. Figure 9 shows that the horizontal tail pro- 
vides more than half of the pitch damping over all but a small range of angle of attack 
between 30° and 33°. In this angle -of -attack range the horizontal tail is probably in the 
low dynamic pressure stalled wake of the wing. At the larger amplitude of oscillation 
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this effect is less pronounced as would be expected because at either end of the oscilla- 
tion the horizontal tail leaves the stalled wake effects of the wing. 

There appear to be no appreciable frequency effects on Cmq + up to an angle 

of attack of 28° as shown in figure 7; increasing the angle of attack beyond approximately 
28° as shown in figure 9 resulted in no consistent trend in the frequency effects on this 
out-of-phase derivative. For the in-phase derivative, CmQ, - k2Cm^ (fig. 8), increasing 
the frequency generally caused this derivative to become more stable, and these effects 
were most pronounced in the angle- of -attack range where static instability had been shown 
to occur. This result is probably related to the well-known fact that flow separation is 
delayed by increased frequency. 

The effects of increasing pitch amplitude are mainly smoothing effects, as would be 
expected, because the nonlinearities in the data would tend to be averaged out with larger 
amplitude motions. Thus, the sharp peaks in both the out-of-phase and in-phase data at 
low amplitudes are flattened out as amplitude is increased (figs. 7, 8, and 9). No consis- 
tent trends were observed in the out-of-phase normal-force derivatives at all amplitudes 
and frequencies, and, consequently, no positive comments as to the effects of amplitude or 
frequency can be offered. (See figs. 7 and 9.) The in-phase normal-force derivatives 
(fig. 8) were very consistent in that the effects of increasing frequency were to increase 
CNq! " k2C]sf^ in the separated flow region. Essentially no effects of frequency occurred 
on either the in-phase or out-of-phase axial-force derivatives. Since there was no appre- 
ciable effect of frequency on either the out-of-phase or in-phase pitching derivatives, the 
data in figure 9 are faired only for the midfrequency range in order to present a clear 
picture of the contribution of the horizontal tail. 

Rolling .- The effects of frequency and amplitude on the out-of-phase and in-phase 
rolling derivatives over the angle -of- attack range of -8° to 112° are presented in fig- 
ures 10 and 11, respectively. 

The damping in roll + Ci^ sin o?j decreases significantly in the stall range and 
at the low oscillation frequencies becomes unstable. This marked reduction of damping 
in roll at the stall is believed to be the basic factor producing lightly damped lateral oscil- 
lations or "wing rock" at high angles of attack as discussed in reference 4. As can be 
seen, the magnitude and sense of damping in roll is a strong function of both amplitude 
and frequency near the stall. As pointed out in reference 4, reduced damping in roll 
will result in a lightly damped Dutch roll motion. Extremely large positive values of 
Cnp + sin a (fig. 10) were experienced in the flat spin range {a = 80° to 90°). Nei- 
ther the source nor significance of these large moments is known. 

Yawing .- The out-of-phase yawing derivatives for three frequencies and two ampli- 
tudes for angles of attack of -12° to 108° are presented in figure 12 for the complete 
model. These data indicate that neither frequency nor amplitude had any significant 
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effect in the unstalled angle -of- attack range, but at higher angles of attack (a 
greater than about 20°) there were considerable effects of these factors, especially in 
Cnr " cos a and Cy^ - Cy^ cos a . The model has stable damping in yaw 

values of Cn^. - Cn^ cos aj through most of the angle- of- attack range, but at spinning 
angles of attack of 60° and greater the model has unstable or positive values of damping 
in yaw. These positive values are an indication that propelling or autorotational yawing 
moments are present on the model, and their effect on the flat spin is discussed in ref- 
erence 5. Figure 13 shows the contributions of the horizontal and vertical tails to both 
the out-of-phase and in-phase yawing derivatives. The vertical tail begins to lose its 
effectiveness at an angle of attack of about 25°. The vertical- and horizontal- tail combi- 
nation becomes effective again, in an adverse sense, in the flat spin range {a between 
70° and 90°) and produces large propelling moments. Removing only the horizontal tail 
completely eliminates the unstable condition and removing only the vertical tail eliminates 
the condition at angles of attack above 75°. The unstable condition is caused by aerody- 
namic interference between the horizontal tail and the vertical tail. This phenomenon is 
explained and illustrated with smoke flow studies in reference 5. Inasmuch as the pro- 
pelling moments were known to be caused by interference between the vertical and hori- 
zontal tail and previous studies had shown large effects of horizontal-tail deflection angle 
on these moments, additional tests were conducted to determine the effect of 6h on 
Cnr “ Cn^ cos a. The results of these tests are presented in figure 14. In the flat spin 
range of 80° to 90° angle of attack the full-up elevator deflection changed Cnr “ cos a 
from propelling to damping, although it is less effective than removing the horizontal tails 
completely. At angles of attack above 95° this combination produced unstable damping 
in yaw. 



CONCLUDING REMARKS 

As a result of the analysis of the oscillatory dynamic stability derivatives of the 
model over an angle -of -attack range of -10° to 110° and at a Reynolds number of 
1.33 X 10^, the following significant points are noted; 

1. The model exhibited stable damping in pitch over the angle-of-attack range 
investigated, but marked reductions in damping in roll were measured at the stall and 
unstable values of damping in yaw were measured at angles of attack associated with 
flat spins. 

2. The model had unstable values of damping in yaw. These values constitute pro- 
pelling moments in a spin and were produced by aerodynamic interference between the 


11 



vertical and horizontal tails. Either removal of the horizontal or vertical tail or full up 
deflection of the horizontal tail eliminated the propelling moments. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., September 25, 1970. 
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TABLE I.- DIMENSIONAL CHARACTERISTICS OF AIRPLANE 


Overall length 

Wing: 

Span 

Area (including leading-edge extension) . . . 

Root chord 

Tip chord 

Mean aerodynamic chord, c 

Leading edge of c rearward of leading edge 

of root chord 

Aspect ratio 

Taper ratio 

Sweepback of 25-percent chord 

Dihedral (inboard 69.5 percent b/2) 

Dihedral (outboard 69.5 percent b/2) . . . . 

Incidence 

Airfoil section: 

Root 

Tip 

Aileron: 

Area (one side) rearward of hinge line . . 
Span (one aileron) (from 44.5 to 

67.0 percent b/2) 

Inboard end chord (base line 

103.24 in. (262.23 cm)) 

Outboard end chord (base line 

155.44 in. (394.82 cm)) 

Spoilers: 

Area (one side) 

Span (from 45.3 to 67.0 percent b/2) . . . 

Inboard end chord 

Outboard end chord 

Horizontal tail: 

Area (in chord plane) 

Movable area 

Span 

Mean aerodynamic chord of horizontal tail . 
Aspect ratio 


57.59 ft (17.55 m) 

38.41ft (11.71m) 

538.34 ft2 (50.01 m2) 

282.00 in. (716.28 cm) 

47.00 in. (119.38 cm) 

192.50 in. (488.95 cm) 

110.76 in. (281.33 cm) 

2.82 

0.167 

45.00° 

0° 

12 . 00 ° 

1 . 00 ° 

NACA 0006.4-64 (modified) 

NACA 0003.0-64 (modified) 

13.08 ft2 (1.22 m2) 

. 4.35 ft (1.33 m) (22.5 percent b/2) 

37.81 in. (96.04 cm) (21.3 percent c) 

34.38 in. (87.33 cm) (27.6 percent C) 


5.44 ft2 (0.50 m2) 
. 4.19 ft (1.28 m) 
. 1.39 ft (0.42 m) 
. 1.04 ft (0.32 m) 

94.70 ft2 (8.80 m2) 
77.40 r 2 (7.19 m2) 
17.705 ft (5.40 m) 
6.143 ft (1.87 m) 
3.30 
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TABLE I.- DIMENSIONAL CHARACTERISTICS OF AIRPLANE - Concluded 


Taper ratio 

Sweepback of 25-percent chord 

Dihedral 

Root chord (at airplane center line) 

Tip chord (theoretical) 

Airfoil section: 

Root (airplane center line) 

Tip (theoretical) 

Hinge-line location, percent Cj^ 

Vertical tail: 

Area 

Span 

Taper ratio 

Root chord 

Tip chord 

Sweepback of 25-percent chord 

Airfoil section: 

Root 

Tip 

Rudder: 

Area (rearward of hinge line) 

Hinge-line location, percent of water-line chords 


0.20 

35.50° 

-23.00° 

. . .107.00 in. (271.78 cm) 
. . . . 21.40 in. (54.36 cm) 

NACA 0003.7-64 (modified) 
NACA 0003.0-64 (modified) 
41.00 

. . . 67.50 ft2 (6.27 m2) 

6.38 ft (1.94 m) 

0.227 

. . 207.15 in. (526.16 cm) 

. . 47.10 in. (119.63 cm) 
58.30° 

NACA 0004.0-64 (modified) 
NACA 0002.5-64 (modified) 

.... 11.07 ft2 (1.03 m2) 
80.00 
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Figure 3.- Photograph of the 0.13-scale model of the twin-jet fighter airplane. 
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(b) Rolling setup. 

Figure 4.- Sketches of test setups for oscillation force tests. 







Q] 


Horizontal 
tails ' 

k 

: - 0 0. 046 

Horizontal q:;";;:— - 
tails 


k 

Horizontal 

tails 

On .. 

□ 

.065 

On 

— □ .065 

On 

On 

On 

A 

.065 

Off --q-lEEETF: 

— <0 


On 

Off 


(a) A0= ± 2 . 470 . 


(b) A0= ±4.910. 


(c) A0= ±8.160. 


Figure 7.- Effect of amplitude and frequency and effect of horizontal tails (at one value of frequency only) 

on the out-of-phase pitching derivatives. Sn = 0*^- 
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because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs^ 
Technology Utilization Reports and 
Technology Surveys. 
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SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washinston, D.C. 20546 


